Distilling Gaussian states with Gaussian operations is impossible 
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We will show that no distillation protocol for Gaussian quantum states exists that relies on (i) arbitrary local 
unitary operations that preserve the Gaussian character of the state and (ii) homodyne detection together with 
classical communication and postprocessing by means of local Gaussian unitary operations on two symmetric 
identically prepared copies. This analysis shows that unlike the finite-dimensional case, where entanglement 
can be distilled in an iterative protocol using two copies at a time, there is no such procedure in the case of 
continuous variables for Gaussian initial states and the above Gaussian operations. The ramifications for the 
distribution of Gaussian states over large distances will be outlined. We will also comment on the generality of 
the approach and sketch the most general form of a Gaussian local operation with classical communication in a 
bi-partite setting. 



In most practical implementations of information process- 
ing devices sophisticated methods are necessary in order to 
preserve the coherence of the involved quantum states. Even 
the mere preparation of an entangled state of spatially dis- 
tributed quantum systems requires such techniques: once pre- 
pared locally and then distributed, an entangled state will to 
some extent deteriorate from a highly entangled state to a 
less correlated state through the process of decoherence. This 
process can quite obviously not be entirely avoided. How- 
ever, one may prepare and distribute several identical entan- 
gled states, and then apply appropriate partly measuring lo- 
cal quantum operations and classical communication to obtain 
states that are similar to the highly entangled original state. 
This is only possible at the expense that one has fewer identi- 
cally prepared systems or copies at hand, but this is a small 
price to pay. Appropriately indeed, this process has been 
given the name distillation [1], as fewer more highly entan- 
gled states are 'distilled' from a supply of many less entan- 
gled states. It was one of the major successes of the field of 
quantum information science to realize that for two-level sys- 
tems such a distillation procedure may be performed on only 
two copies at a time, and it requires only two steps: (i) a local 
unitary operation, and (ii) a local measurement, together with 
the classical communication about the measurement outcome. 
Based on the measurement outcome further local unitary op- 
erations are then implemented. 

Such distillation protocols may also be of crucial impor- 
tance in the infinite-dimensional setting. Quantum informa- 
tion science over continuous variables has seen an enormous 
progress recently, both in theory and experiment, mostly in- 
volving Gaussian states of field modes in a quantum optical 
setting [2-A]. Quite naturally, one should expect that a simi- 
lar distillation procedure also works for Gaussian states in the 
infinite-dimensional case, also under the preservation of the 
Gaussian character of the state. If one transmits two pure two- 
mode squeezed Gaussian states through lossy optical systems 
such as fibers, the corresponding modes being from now on 
labeled Al, A2, Bl, and B2, one obtains two identical copies 
of less entangled symmetric states [5]. A feasible distillation 
protocol preserving the Gaussian character may consist of the 
subsequent steps (see Fig. 1): 



(i) Application of any local Gaussian unitary operation. 
That is, one may implement any unitary operations Ua and 
Ub on both Al and A2 on one hand and Bl and B2 on 
the other hand corresponding to symplectic transformations 
[6] S A ,S B G 5p(4,R) [7]. This set includes all two-mode 
and one-mode squeezings, mixing at beam splitters and phase 
shifts. To specify these operations 20 real parameters are nec- 
essary. Note that we do not require both parties to realize the 
same transformation. 

(ii) A homodyne measurement on the modes A2 and B2. 
The parties communicate classically about the outcome of the 
measurement, and may postprocess the states of modes Al 
and Bl with unitary Gaussian operations. 
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Al A2 Bl B2 

FIG. 1. The class of considered feasible distillation protocols. 

The main result of this Letter is that very much as a sur- 
prise, none of these protocols amounts to a distillation pro- 
tocol. No matter how ingeniusely the local unitary operation 
is chosen, the degree of entanglement can not be increased. 
The optimal procedure is simply to do nothing at all, which 
means that at least no entanglement is lost [8]. The degree of 
entanglement will be measured in terms of the log-negativity, 
which is defined as -EW(p) = log 2 \\p TA || for a state p, where 
|.| denotes the trace norm, and p TA is the partial transpose 
of p. The negativity has been shown to be an entanglement 
measure in the sense that it is non-increasing on average un- 
der local operations with classical communication [9], and is 
to date the only known feasible measure of entanglement for 
Gaussian states. For pure (and for symmetric mixed) Gaus- 
sian states it is related to the degree of squeezing in a mono- 
tone way (see, e.g., [10]). This means that as a corollary of 
the main result, it follows that with Gaussian operations as 



1 



specified above one cannot transform two identically prepared 
two-mode squeezed vacua into a single two-mode squeezed 
vacuum state with a higher degree of squeezing. 

We will start by fixing the notation. Gaussian states [11] 
of an n-mode system are completely characterized by their 
first and second moments. The first moments are the expec- 
tation values of the canonical coordinates. The second mo- 
ments can be collected in the real symmetric covariance ma- 
trix T G C(2n) C M(2n,M), where M(2n,E) denotes the 
set of real 2n x 2n-matrices, and C(2n) the subset of matri- 
ces obeying the Heisenberg uncertainty principle [11]. The 
linear transformations from one set of canonical coordinates 
to another which preserve the canonical commutation rela- 
tions form the group of real linear symplectic transformations 
Sp(2n,R) [6]. A symplectic transformation S changes the 
covariance matrix according to T i — ► STS T , while states un- 
dergo a unitary operation p i — ► U pU^ . The n = A modes Al, 
A2, Bl, B2 will be equipped with the canonical operators 
(Xai, Pai, Xb2, Pbi)- To make the notation more trans- 
parent, both tensor products and direct sums will carry a label 
indicating the underlying split, meaning either A, B or 1, 2. 
We state the main result of this Letter in form of a theorem: 

Theorem. - Let p® pbe two identically prepared symmetric 
Gaussian states of two-mode systems consisting of the parts 
Al, A2, Bl, and B2, respectively, each of which having the 
covariance matrix 
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let Sa,Sb G Sp(4, M.) be any symplectic transformations 
with associated unitaries U a and Ub, and let 

p' = (Ua ®a,b U b )(p ®i,2 p)(U A ®a,b U b ) ] . (2) 

Then any state p" that is obtained from p' via a selec- 
tive homodyne measurement on systems A2 and B2 satisfies 
En(p") < En(p), that is, the degree of entanglement can 
only decrease. 

The proof of this statement will turn out to be technically in- 
volved, and while the statement itself is concerned with prac- 
tical quantum optics, the techniques used in the proof will be 
mostly taken from matrix analysis [12]. In order to give the 
general argument more structure, the proof is split into sev- 
eral lemmata. The entire proof will be formulated in terms of 
covariance matrices, rather than in terms of the states. 

The log-negativity of a state a of a two-mode system can 
be easily expressed in terms of the entries of the associated 
covariance matrix 7 G C(4). The latter can be partitioned in 
block form according to 

7 = ( It Is ) ' 7A ' 7s ' 7c e M(2 ' R) - (3) 

The log-negativity En(ct) is then given by [9] 



-(logo/)( 7 )/2, if/( 7 )<l, 
otherwise. 



(4) 



where the function / : C(4) — > M+ is defined as 

/( 7 ) - ((det[ 7 A] + det[ 7B ])/2 - det[ 7 c]) 

- [((detbA] + det[7u])/2 - det[ 7c ]) 2 - det[ 7 ]] 1/2 . (5) 

The covariance matrix associated with the Gaussian state p' 
in the Theorem will be denoted as T' G C(8). For any 
Sa,Sb S Sp{A, R) this covariance matrix of the modes Al, 
A2, Bl, and B2 becomes V := (S A ®a,b S B ){T (a) ©1,2 
r(°))(S' J 4 ®a,b Sb) t The first step is to relate the covariance 
matrix V" associated with the state after the measurement to 
a Schur complement [12]. This Schur complement structure 
is a general feature of Gaussian operations and will be further 
discussed at the end of the letter. 

Lemma 1. - Let V e C (8) be a covariance matrix of systems 
Al, A2, Bl, and B2 associated with a state p', which can be 
written in block form as 



V 



C\ C3 
C2 



(6) 



where Ci,C%,C3 G M(4,R). The covariance matrix of 
the state that is obtained by a projection in A2 and B2 
on the pure Gaussian state with covariance matrix Dd := 
diag(l/d, d,l/d,d) G C(4), d > 0, is then given by 



M d = C 1 -C 3 {C2 + ti 2 d )- 1 C% 



(7) 



Proof. This statement can be most conveniently be shown 
in terms of the characteristic function \ [HL By employing 
the Weyl (displacement) operator, the state p' associated with 
the covariance matrix V can be written in terms of the charac- 
teristic function according to p' = (I/71- 4 ) / d 8 £ W(— 0x(0 
(see, e.g., Ref. [13]). The projection corresponds on the level 
of the characteristic function therefore to an incomplete Gaus- 
sian integration. The characteristic function associated with 
the modes Al and Bl can then be written as 

x(a , ...,&)=/ ^-** e -er>v2 e -&&+&-£&+® 

= \C 2 + D 2 |- 1 /2 e -«i,-.€4)r"( €l ,...,€4) r /2 ; (g) 

with Md defined as in Eq. (7). ■ 
Hence, the resulting covariance matrix is given by the Schur 
complement C\ - C 3 (C 2 + D 2 -) -1 ^ of the matrix 



r' — 



C\ c 3 

C 2 + Dj 



(9) 



with respect to the leading principal submatrix C\ . The addi- 
tional matrix D\ originates from the projection in the modes 
A2 and B2. Note that although this Lemma has been for- 
mulated in terms of the projection on a certain class of pure 
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Gaussian states, it applies to the projection on any pure Gaus- 
sian state in the modes A2 and B2: the projection on any 
other pure Gaussian state can be realized by an appropriate 
choice of the symplectic transformations Sa and Sb- Ideal 
homodyne detections can now be formulated as projections 
on 'infinitely squeezed' pure Gaussian states [13]. The cen- 
tral feature is that the initial first moments do not affect the 
form of the covariance matrix after the measurement. Lemma 
2 gives the form of the resulting covariance matrix in case of a 
homodyne detection in modes A2 and B2. In the limit d —> 
the matrix Dd gives rise to a projection operator, and the in- 
verse becomes a Moore Penrose inverse (MP) [12]: 

Lemma 2. - In the notation of Lemma 1, the covariance ma- 
trix of modes Al and Bl after a selective homodyne measure- 
ment in modes A2 and B2 is given by 



T" := lim r d 

d^O 



d -C 3 (7rC 2 7r) M/J C : 



MP r <T 

3 i 



(10) 



where n — diag(l, 0, 1, 0). 

Equipped with these preparatory considerations, we will 
now turn to the core of the proof. In order to be able to evalu- 
ate the logarithmic negativity according to Eq. (4), one needs 
to know the values of the invariants under local symplectic 
transformations, i.e., the determinants of four submatrices. 
To find an expression for all these determinants is however 
a quite difficult task. Instead, we will later make use of an 
upper bound of the logarithmic negativity that only involves 
determinants of principal submatrices [12] of T". 

Lemma 3. - Let T" G C(4) be defined as in Lemma 2. Then, 
independent of S At Sb £ Sp(4:,~Rl), 



det[T"] = det[Y {0) ] = (a 2 



„2\2 



(ID 



Proof. According to Lemma 2, T" is given by T" — 
linid^o M d . The Schur complement of the matrix T' d as de- 
fined in Eq. (9) is related to T' d and one of its principal sub- 
matrices via the congruence 



U X 

o u 



1 A 



u o 



r„ 



o 



C 2 + D 2 



where X := —C?,(C2 + D 2 l )~ 1 . Hence, according to 
the determinant multiplication theorem we obtain det[r^] = 
det[r d ] det[C" 2 + D% which yields in the limit d -> 

det[r"] = det[Pr^P + (1 8 - P)]/det[Qr^Q + (1 8 - Q)], 

where the projections P and Q are defined as P := 
diag(l, 1,1, 1,0, 1,0,1) and Q := diag(0, 0,0, 0,0, 1,0,1). 
With these tools, it is feasible to directly prove the state- 
ment of Lemma 3 by parameterizing Sa,Sb G Sp(4,~R). 
Every S G Sp(A, R) can be written as a product S = 
VDW, where V, W e Sp(4,R) n 5*0(4), and D := 
diag(di, l/d\, d 2 , l/d 2 ) with di, d 2 G R [6]. ■ 
Lemma 4. - Let T" G C (4) be defined as in Lemma 2, and let 
T' A and T B be the principal submatrices belonging to mode 
Al and Bl. Then, for all S A , S B & Sp(4, R), 



det[T" A ] < det[T { l ] ] = a 2 , det[T B ] < det[T ( ° ] ] = a 2 . (12) 

Proof. T" is defined as the covariance matrix correspond- 
ing to modes Al and Bl after the projective measurements in 
both A2 and B2. Let us assume that one first performs the pro- 
jective measurement in A2, leading to a the covariance matrix 
-/V4 G C(2) of the reduced state of Al. The covariance matrix 
T' A after the projection in B2 is then obtained as a Schur com- 
plement. In particular, T' A can be written as T' A = Na - R, 
where R G M(2, R) is a real symmetric positive matrix. 
Hence, as T' A and Na are also positive, det[T^] < det[7V^] 
[12]. In other words, one obtains an upper bound for det[r^] 
when considering only a projective measurement in Al. The 
statement of Lemma 4 follows from Lemma 3 in the special 
case that c = 0: one can after a few steps conclude that then 
det[iV A ] = a 2 , independent of S A , S B € Sp(A, R). The same 
reasoning applies to T' B . ■ 
The most important step is now an appropriate upper bound 
of the log-negativity of the resulting state. The actual bound 
might appear somewhat arbitrary, but it will turn out that it is 
exactly the tool that we need in the last step of the proof. 

Lemma 5. - Let 7 G C(4), partitioned as in Eq. (3). Then 

/(7) > .9(7) := [((det[ lA ] + det[ lB ]) /2) 1/2 (13) 
- {{det[ lA ]+det[ lB ])/2 - <fef[ 7 ] 1/2 ) 1/2 ] 2 - 

Proof. g(j) can be expressed in terms of / as 5(7) = f {"{'), 



Ja 7c 



1'a 



1'b 



a 1 1 



2, 7c 



d 
-c' 



where a! := ((det[ 7j4 ] 2 + det[ 7B ] 2 )/2) 1 / 2 and c' := (a' 2 - 
det[ 7 ] 1/2 ) 1/2 . Hence, one has to prove that 7(7') < /( 7 ). 
Firstly, note that det[ 7 ] = det[ 7 ']. Secondly, (det[7^] + 
det[ 7 s])/2 = a' . Therefore, it remains to be shown 
that c' 2 > |det[ 7 c]|. This inequality is equivalent with 
[(det[ 7A ] +det[7 B ])/2 - |det[ 7c ] |] 2 -det[ 7 ] > 0, which is a 
valid inequality, as 7 G C(4). ■ 
Proof of the Theorem. Let V" G C(4) be the matrix de- 
fined as in Lemma 2. The log-negativity of the correspond- 
ing state of modes Al and Bl is given by — (log of)(T"), 
if the final state is entangled at all, as we will assume from 
now on. Lemma 5 yields the bound f(T") > g(T"). In 
<?(r"), however, only the determinants of the principal sub- 
matrices are needed, bounds of which are available by virtue 
of Lemma 3 and 4. The function h : [y, 00) — > R + with 
h(x) = {x 1 / 2 — [x — y) 1 / 2 ) 2 , y > 0, is a strictly mono- 
tone decreasing function of x. Therefore, using Lemma 3 
and 4 one can conclude that <?(r") > g(T^). Moreover, 
ff (r(°)) = /(r(°)), due to the special form of r(°), as can be 
easily verified. Hence, f(T") > g(T^) = f(T^), which 
leads to -(logo/) (r") < -(logo/)(r(°)). This is finally 
the desired result: it means that the degree of entanglement 
can only decrease. ■ 
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We will finally comment on the generality of the approach. 
A general Gaussian operation is a quantum operation that 
maps all Gaussian states on Gaussian states [3]. Any gen- 
eral Gaussian local operation with classical communication 
(LOCCG) - trace-preserving or non-trace-preserving - can be 
decomposed into the subsequent steps: (i) Appending locally 
additional modes that have been prepared in a Gaussian state 
[3]. (ii) Application of any local unitary Gaussian operation 
on both the original and the additional system. These com- 
prise operations corresponding to symplectic transformations 
and displacements in phase space, (iii) Projections on pure 
Gaussian states or ideal homodyne detections, which give rise 
to Schur complements on the level of covariance matrices as 
described above, together with the classical communication 
about the outcome (real numbers in case of homodyne detec- 
tion, bits in case of dichotomic measurements including the 
projection on a pure Gaussian state in one outcome), (iv) mix- 
ing, such that the resulting state is Gaussian, and (v) a par- 
tial trace, which corresponds to considering a certain princi- 
pal submatrix of the covariance matrix only [14]. The proof 
is therefore restrictive in the sense that only two copies at a 
time are considered, other projections on Gaussian states are 
excluded, and no additional modes are allowed for. The state- 
ment of the present paper proves that iterative protocols in 
strict analogy to the corresponding methods in finite dimen- 
sional settings certainly do not work. Indeed, the findings 
strongly suggest that Gaussian states cannot be distilled at 
all with Gaussian operations. Then (less feasible) non-linear 
physical effects [15] would have to be made use of in order 
to distill from a supply of Gaussian two-mode states [16]. 
Such techniques would then also be necessary for the realis- 
tic implementation of quantum repeaters [17] for continuous- 
variable systems when it comes to the distribution of highly 
entangled Gaussian states over large distances. 
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